Muon Dynamics in a Toroidal Sector 

Magnet 

Juan C. Gallardo 1 , Richard C. Fernow and Robert B. Palmer 

Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 



Abstract. We present a Hamiltonian formulation of muon dynamics in toroidal sector 
solenoids (bent solenoid). 



I INTRODUCTION 



The present scenario for the cooling channel in a high brightness muon collider 
[1] calls for a quasi-continuous solenoidal focusing channel. The beam line consists 
of a periodic array of hydrogen absorbers immersed in a solenoid with alternating 
focusing field and rf linacs at the zero field points. 

The simple ^ energy loss in conjunction with multiple scattering and energy 
straggling leads to a decrease of the normalized transverse emittance. Reduction 
of the longitudinal emittance could be achieved by wedges of material located in 
dispersive regions; at least in principle, this scenario seems appropriate to obtain 
effective 6-D phase space reduction. [2] 

A conventional chicane is a dispersion element but its use presents a serious 
challenge, as it is very difficult to integrate it with the solenoidal channel. The 
matching into the periodic solenoidal system imposes constraints on the Twiss 
parameters of the beam which seems not easily achievable. A possible alternative 
is the use of curved solenoids in conjunction with wedge absorbers as suggested by 
one of the authors. [3] 

Solenoids and toroidal sectors have a natural place in muon collider design given 
the large emittance of the beam and consequently, the large transverse momentum 
of the initial pion beam or the decay muon beam. Bent solenoids as shown in Fig.l 
were studied for use at the front end of the machine, as part of the capture channel 
[4] and more recently as part of a diagnostic setup to measure the position and 
momentum of muons. [5] 
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II TOROIDAL SECTOR SOLENOID 



If we restrict ourselves for the moment to a horizontal bending plane, the mag- 
netic field inside of the solenoid and near the axis has a gradient (field lines are 
denser at smaller radius) described approximately by B(x,y,s) ~ B s e s with 

B s {x,y,s) ^ 



(1 + hx) 

where s is the coordinate along the particle trajectory and h — jj- is the curvature 
at the position s, with R Q the radius of curvature. As a consequence of the curvature 
of the trajectory and the corresponding magnetic gradient, the center of the particle 
guide orbit, averaged over the Larmor period, drifts in a direction perpendicular to 
the plane of bending [6]. The combined drift velocity can be written as, 

dr B rrif, 2 2 (R Q x B) 

Tt =V *B + ^ i2V * +V±) (l + hx)E>B* - (2) 
and the magnitude of the transverse drift velocity is 

Clearly a y-position versus energy (v\\) correlation will develop as the muon beam 
travels along the toroidal sector solenoid. 

From Eq.2 above we notice that if we include an additional vertical field, a dipole 
with a curvature equal to that of the bent solenoid for the reference energy, i.e. 

Bd ~ -— vjfcjftfy ( 4 ) 



then Eq.2 reduces to, 



dr B s _ , , n B 



and consequently, particles with the chosen energy will not drift vertically and will 
remain on the axis of the bent solenoid. Those particles with larger energy will 
drift upward (positive y-direction) and those with lower energy downward (negative 
y-direction), achieving the needed dispersion. The magnitude of the dispersion is 
given by 



D y = = 2«h(^\ (6) 

y q Bl \qB s ' 



where p is the chosen momentum corresponding to zero dispersion. 
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FIGURE 1. Schematic of emittance exchange system with a bent solenoid 



III DYNAMICS IN A TOROIDAL SOLENOID. 
HAMILTONIAN FORMULATION 

From general principles of classical mechanics and following the usual approxi- 
mations in accelerator physics [7] the normalized Hamiltonian reads, 

H(x,p x ,y,p y ,z,S; s) ps -^(1 + hx)A s (s) - (1 + 8)hx 

+ SttS{^-^) 2 + ^-^) 2 } (7) 
where the path length s is the independent variable, p x , p y are normalized momenta 



with respect to po, the initial reference total momentum p = ^Jp 2 . + p 2 , + P 2 ', z — 

s — P ct, 5 = and A = (A x ,A y ,A s ) is the vector potential. The vector 

potential satisfies the gauge condition V • A = 0. 

In the accelerator frame of reference, i.e. the Frenet-Serret coordinate system 
defined by the metric, 

da 2 = dx 2 + dy 2 + (1 + hxfds 2 (8) 
the equations for the coordinate unit vectors are 

^ = h(s)e s , ^ = , ^ = -h(s)e x (9) 
ds ds ds 



The magnetic and electric fields are obtained from 



B = V x A and E = c(3 Q 



dA 

dz 



(10) 



with 



V x A 



(i+hx) hx)A s — d s Ay\ e x 

+ JT+Kx) i d sA x - d x (l + hx)A s } e y 
+ \d x Ay dyA x } e s 



and 



V • A 



1 + hx) 















— (A x (l + hx)) + — (A y (l + hx)) + —A 



dx 



dy 



ds 



'ir 



(12) 



The lowest order approximation for a toroidal solenoidal field is given by Eq.l. 
The corresponding vector potential in the next order is [8], 



A 



-Ib. 



y 



B n 



2^(l + hxf x+ 2K ln{1 + hx) ^ 



and the corresponding magnetic fields are: 



B x 



i 



i 



(l+h(s)x) 2h(s) 



{^\n(l + h(s] 



x 



_Bo_h' 

h(s) 



h'(s)]*(l + h(s)x) + h>(s) T ^_} 



(13) 



(14a) 



2(1 + h(s)x) 



dB B Q (s)xy 
-y-h (s) 



ds 



1 + h(s)x 



(14b) 
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1 + h(s)x 

One possible second order approximation for the vector potential of a dipole is 



(14c) 



A x 

Ay 

A a 



B D h"{s)x{y 2 -\x 2 ) 
2h(s) (l + h{s)x) 
B D h'(s) xy 
2h(s) (1 + h(s)x) 

-^){ 1 + k{s)x - h ' {s) ^ 2 -^} 



(15a) 
(15b) 
(15c) 



Substituting the total vector potential into the Hamiltonian, and dropping some 
constants we can write, 



H to f. S oi.( x ^Px,y,Py,^,S; s) « \(hx) 2 -bhx 
+ [ Ww{(p X -^ x f + (p y -^A y r} (16) 

We have written a simple Fortran program to solve the equations of motion from 
the above Hamiltonian ; its result for a few representative cases of interest are 
shown in Fig. 2. 




A second order expansion of the bent solenoid magnetic field given in Eqs. 14 has 
been used together with a second order expansion of the dipole magnetic field in the 
cooling simulation program ICOOL [9]. Fig.3a shows an example of the dispersion 
D y in a bent solenoid obtained in ICOOL as a function of the dipole strength Bp. 
It is apparent that the dependence of Eq. 6 on Bp is well satisfied. Likewise, Fig 
3b shows simulation results for the dispersion as a function of B~ 2 . Again we see 
that the mentioned equation gives a good representation of the results. 
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